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We consider the six dimensional brane world model, where the brane is described by a localized 
solution to the baby-Skyrme model extending in the extradimensions. The branes have a cosmologi- 
cal constant modeled by inflating four dimensional slices and we further consider a bulk cosmological 
constant. We focus on the topological number three solutions and discuss the localization mecha- 
, nism of the fermions on the above 3-branes. We discuss interpretation of the model in term of quark 

. third generation mass and in terms of the inflation history. We argue that the model can describe 

■ various epochs of the early universe by suitably choosing the parameters. We further discuss the 

localization properties of gauge fields on the brane and argue that this is achieved only for specific 
values of the electromagnetic coupling, providing a quantization to the electric charge. 



PACS numbers: ll.lO.Kk, 11.27.+d, 11.25.Mj, 12.39.Dc 
^ , Keywords: Field Theories in Higher Dimensions; Solitons Monopoles and Instantons 

I. INTRODUCTION 

m 
> 

■ Theories with extradimensions have been expected to solve the hierarchy problem and cosmological constant prob- 
, lem. Experimentally unobserved extradimensions indicate that the standard model particles and forces are confined 
f — ■ to a three brane [Jill. The Randall- Sundram (RS) brane model in five space-time dimensions 0,13] shows that the 
04 \ exponential warp factor in the metric can generate a large hierarchy of scales. The brane theories in six dimensions 
. in models of topological object show a very distinct feature towards the fine-tuning and negative tension brane prob- 
lems. In the context abelian strings 0-111 were investigated, showing that they can realize localization of gravity for 
negative cosmological constant. For the magnetic monopoles, similar compactification was achieved for both positive 
and negative cosmological constant [l^ . 

There are two main contexts in which solitons appear in field theories: one is related to the strings and the magnetic 
monopoles in non-abelian gauge theories, and the others are kinds of non-linear type models like the skyrmions, 
hopfions [ill. IT^ . The Skyrme model is known to possess soliton solutions called baby-skyrmions in two dimensional 
space Il3l4l5ll . The warped compactification of the two dimensional extra space by such baby skyrmions was already 
studied [l6[ for negative bulk cosmological constant, based on the assumption that the cosmological constant inside 
the three branes is tentatively set to be zero. Addressing the non-zero cosmological constant inside the branes has 



been considered for case of the strings [17[ and the monopoles [18|. Along these directions, we also have studied the 
baby-skyrmion brane with both positive brane cosmological constant and a bulk cosmological constant [l9| . In this 
paper, we employ these solutions as backgrounds for the fermions. 

Study of fermion and gau ge fields localization on topological defects have been extensively studied with co-dimension 
one [2Cl| - l27j and two [284371. M any years ago, particle localization on a domain wall in higher dimensional space time 
was already addressed 3^, [S^] . The authors suggested the possibility of localized zero- modes of fermions on the one 
dimensional kink background in 4-1-1 space-time with Yukawa-type coupling. Later, localization of chiral fermions 
on RS scenario was discussed in [20j . Analysis for the massive fermionic modes was done by Ringeval et.al., in 
[2l| and later several studies have followed [25l - [27j . For co-dimension two, the localization on higher dimensional 
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generalizations of the RS model was studied within the coupling of real scalar fields [28| . Many studies have followed 
and most of them arc based on the Abelian Higgs, or Higgs mediated models with the chiral fermions. 

The problem of mass hierarchy in the Standard Model (SM) fermions [13] has been discussed in many articles 
based on the brane worlds [23. [27l - [30l [32l. [ssl [sgI IstI. |41| - |4^ in several mechanisms. For example, in [s^ the fermions 
have quantum numbers of the rotational momenta which are origin of the generation of fermions. The authors of 
[3^ [ssj deal with this problem with somewhat different approach. Conical singularity of the background branes and 
the orbital angular momentum of the fermions around the branes are the key role for the generation. In [29I |43 | 
hierarchy between the fermionic generations are explained in terms of multi- winding number solutions of the complex 
scalar fields. They observed three chiral fermionic zero modes on a topological defect with winding number three and 
finite masses appear the mixing of these zero modes. Although any discussion of brane construction is absent in their 
discussion, the idea is promising. In [4^, the authors have taken into account more realistic standard model charges. 

Among these, we consider the localization of the fermions on the "inflating " baby skyrmion branes with topological 
charge three. The localized modes of fermions are confirmed through the analysis of spectral flow of the one particle 
state [i^ . According to the Index theorem a nonzero topological charge implies the zero modes of the Dirac operator 
[46| . The zero crossing modes are found to be the localized fermions on the brane. So the generation of the fermions 
is deflned in terms of the topological charge of the skyrmions with a special quantum number called grandspin K^^. 
There arc different proflles of the zero crossing behavior for different K^, and it is the origin of the flnite mass in 
our point of view. Our main concern of this paper is localization of the fermions inside the brane with positive 
cosmological constant, i.e. the inflating branes. 

We further try to make more sense of such models from a phenomenological point of view. Indeed a four dimensional 
cosmological constant is now admitted, and the model under investigation has positive curvature. However, the 
positive curvature is here modeled by inflation, which is clearly suitable for early universe bit might be inaccurate 
for present epoch. We argue that suitable choice of parameters might still give a reasonable approximation of the 
underlying physics thanks to the fact that only he Hubble parameter enters the (background) equations. Further 
changes affect the details of the four dimensional internal physics without spoiling the localization properties. The 
present model can be suited to describe more general positive curvature four dimensional branes. The Hubble 
parameter was quite large value at the inflating period; it dramatically reduced and is almost negligible at the present 
epoch. Interestingly, the history of such change of the Hubble parameter makes a new spectral flow behavior which 
may explain the particle creation in our Universe. 

Next, as a straightforward extension of the study, we consider localization of gauge fleld minimally coupled to the 
braneworld. It is certainly natural to take into account of particle on the brane after the fermions. As a result, we 
could treat the quarks (and/or electrons) and also the photon on the brane which certainly ensures to get closer to 
the Standard Model contents. We discuss the localization properties of the gauge fleld as well as its effective four 
dimensional interpretation. Our flndings indicate that the four dimensional fleld is localized only for very speciflc 
values of the electromagnetic coupling. 

Note that we do not consider the backreactions of the brane due to the gauge flelds nor the fermions. Somehow, we 
assume that the matter content of the brane has a negligible contribution to the brane energy, i.e. that the content 
is low energy and the extradimensional physics is of high energy, which seems reasonable. 

This paper is organized as follows. In the next section we briefly describe the Einstein-Skyrme system in six 
dimensions. Several types of solutions are shown in Sec. II. We discuss with the formalism of the fermions coupling to 
the skyrmion solutions in Sec. HI. The spectral flow analysis is introduced in this section. In Sec IV, we consider the 
effects of the inflation to the spectra of the fermions. In Sec. V, we discuss the localization of massless electromagnetic 
flelds in the background of the brane. We summarizes in Sec. VI. 



Let us flrst introduce the notations and remind the model we study. The bosonic sector is the same as in our 
previous study so we briefly summarize it. The total action for the gravitating baby-Skyrme model is of the form 
S := S'grav + iS'baby The gravitational part of the total action 



II. THE GRAVITATING BABY-SKYRME MODEL IN SIX DIMENSIONS 



A. 



The model 




(2.1) 
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is the generalized Einstein-Hilbcrt gravity action, where A(g) is the bulk cosmological constant, and X(6) = 87rG(6) 
87r/M^^g-j. The second term S'baby in the total action stands for the baby-Skyrme action and is given by 



^1 



baby 



I d'^xy^[^{dMn) ■ {d''n) ~ ^{dMU X dNu f ~ KoV{n)] . (2.2) 



Here n is a scalar triplet subject to the nonlinear constraint n ■ n = 1, and V{n) is the potential term with no 
derivatives of n. The coefficients ^0,2.4 in Eq. p.2p are the coupling constants of the gravitating baby-Skyrme model. 
Assuming axial symmetry for the extradimcnsions, the metric can be written in the following form 

ds^ = Af2(p)gWrf^Mrfa.- _ _ l^(p)de^ (2.3) 

where p G [0, 00) and 6 G [0, 2tt] are the coordinates associated with the extra dimensions. 

We further model a cosmological constant on the brane by considering the following form of the four dimensional 
subspace (described by g^^J in Ea. (|2.3p ) 

ds^4) = gf^dx^'dx" = dt^ - S,je^"°^dx'dx^ (2.4) 

where the constant Hq is so called the Hubble parameter. 

The ansatz for the scalar triplet n is given by the hedgehog Ansatz: 

n = (sin f{p) cos(n 6*), sin f{p) sm{n 9), cos f{p)). (2-5) 

Let us note that there are some variations [13] for choosing the potential term V{n) in the baby-Skyrme model 
(|2.2|) . Here we use the so-called old baby skyrmions potential, which reads 

y(n) = l-n-n(°°) = l + cos/(p), (2.6) 
where n^°°'' = limp_>oo '>^{Pi is the vacimm configuration of the baby-Skyrme model. 



B. Field equations, boundary conditions and parameters of the model 

The equations of motions are extensively described in [l^ so we skip the details. Here we only present the definitions 
of the reduced variables for the sake of the readers' understanding: 

^p, L{r):^,[^l{r\ (2.7) 
K4 V K4 



and dimensionless parameters 



K4 K4 K4 2 

a ■■= X(6)K2, P := A(6) — , M := kq — ,7 = —Hq 

^2 ^2 



Note that it is possible to interpret 7 as a positive cosmological constant in the four dimensional subspace of the full 
model using the four dimensional effective theory following the lines of Indeed in this case, g^^J is such that 

G^^J = ZH^g^^fl , where G^^^ is the Einstein tensor computed with g^'^-'. Note that replacing the four dimensional 

subspace by another Einstein spacetime satisfying g'"^^ = SH^g'^J, such as the Schwarzschild-de Sitter spacetime, 
leads to the same equations. 

Let us remind the boundary conditions too: 

/(0) = -(m-l)^, /(oo) = ^, (2.9) 

where m, n S Z, for the baby Skyrme field and 

L(0) = 0, L'(0) = = 1, M(0) = 1, M'(0) = = 0, (2.10) 

dr r=o dr 



r=0 



for the metric fields. The above boundary conditions are required for regularity and finitcness of the energy. 

Remember that considering the hedgehog Ansatz (|2.5|) with the boundary condition (j2.9|) leads to a topological 
charge (or winding number) given by 



N 



If n 

— / n • {dpu X dgn)dpde + (-1)^"'"] e Z. (2.11) 



Another useful quantity is the rescaled Ricci scalar which will be used later and is given by 

2L" , WM' , MI" , l2{M'f 127 
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(Q-,/3,7,;j,m,n) = (0.005,0.,0.175,0.1,1.,3.) 




FIG. 1: A typical solution for n = 3,m — 1. 



C. Asymptotic solutions 



Here we briefly remind the asymptotic solutions available in the gravitating baby Skyrme model. 
For the large asymptotics r — oo, we can set /(r) = 7r,i.e., the baby-Skyrme field is a topological vacuum. For 
/? = we have a cigar-type set of solutions given by 



Af,(r) =7(r-ro), L^ir) = Lo = -/ Cl 



Periodic solutions are found for /? > 0: 



'IO72 



/a/3 
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^™ V "i7r(^ ^ '^o) , Lp{r) = io cos \ — (r - tq) 
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(2.13) 



(2.14) 



For /3 < we have diverging solutions 



'loy^. 



-a/3 



sinh 



To" 



(r - ro) , Ld{r) = Lq cosh 



-a/3 

To" 



(^^ - ^-o) 



(2.15) 



The near origin behaviour of the functions /, L, M subject to the boundary conditions (|2.9p and (|2.10l) is given by 

^3 



f{r) = -(m-l)7r + /(")(0)r"+O(r)"+\ i(r) = r + Ai- + 0(r)^ 
M[r) =1 +(Ai +7)— + 0(r) , Ai = , 



(2.16) 



where /^-"^ stands for derivative of /. Note that higher order corrections are straightforward to compute. 



D. The = 3 solutions 



According to (|2.1ip . the case of topological number three can be achieved only for n = 3 and odd values of m. 
The simplest case is of course (m, n) — (1, 3) but we also consider some solutions for (m, n) — (3, 3). In the following 
analysis, we mainly focus on the case /3 = 0, i.e. the case of no bulk cosmological constant. 

We numerically solve the system of ordinary differential equations with the solver Colsys [i^. Here we report 
some notable features of our solutions when varying the model parameters. The solutions we built share the common 
qualitative properties as their n = 2 counterpart, up to some details which does not affect the underlying physics. 
For example, there is a maximum value of 7 that becomes a critical point where the solutions cease to exist for each 
/i, a. The baby skyrmion somehow shrinks to smaller radii as 7 grows until it disappears. 

We show the solutions with m = 1 on Fig. [1] The metric function M develops a minimum for certain range of 
parameter a for {m,n) = (1,3). The interpretation is the following: the baby Skyrme wants to shrink the four- 
dimensional slices up to some radii where it exercises its gravitational interaction, but then for larger radius, the 
spacetime is no longer influenced by the brane and the four dimensional slice grows. Interestingly, such phenomenon 
do not occur for m = 3. This is illustrated on Fig. [5] 
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FIG. 2: For some range of parameters, the warp factor M develops a minimum, it is tempting to assume that the solution ceases 
to exist when the warp factor becomes zero. The parameters of the above solutions are 7 = 0.001, /3 = 0, /i = 0.1, m = 1, n = 3 
and the range of a is given on the figure. 

The matter distribution in the case (m,n) = (1,3) is of a disk shape with more matter at the outside. As 7 
increases, even if the root mean square radius decreases, the matter distribution becomes dominant on a ring close to 
the origin but remaining of a disk shape. 

Instead, for m = 3, the derivative of the function / becomes peaky close to the origin as 7 increases (again for fixed 
a, /i. As a consequence, the matter distribution becomes also peaky on a ring close to the origin. Note that in the 
case m = 3, the baby skyrmion has indeed a ring shaped matter distribution and no matter is present at the center. 
This is illustrated on Figs. [3]and|4l 



m=l, 11=3, ci-=.005. 13=0, n=.l, y=.00l m=l, n=3, a=.005, 13=0, ii=.l, y=.030 




FIG. 3: Matter distribution for {m,n) = (1,3) in two cases, the center is not vacuum and the ring shape is more pronounced 
as 7 increases. 



i=3,n=3,a=.001,/?=0,,j=.01,7=.0001 m=3. n=3. a-=.001, /i=0. /i=.Ql. y=.0012 m=3, n=3. a-=.001. /i=0. /i=.Ql, y=.0153 
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FIG. 4: Matter distribution for (m, n) — (3, 3) in two cases, here the center is vacuum and the distribution is of a ring shape. 
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III. MASSIVE FERMIONS 



The fermions coupled with the baby-Skyrnic field have an isospin doublet structure. Furthermore, we tentatively 
assume that the fermions within an iso-doublet are degenerate. Of course, the assumption is not valid especially for 
the heavier quark sectors. In order to recover it, we would have to introduce additional terms which explicitly break 
the symmetry. Or, different choice of the parameters of the baby-skyrmions for each flavor can exert similar effect, 
which wc shall apply. Also, wc ignore inter-gcneration mixing of SM particles. 



A. Basic formalism 



We couple the fermions to the baby Skyrmc background by virtue of the following action 

^fcrmion = J d"^ x^imT"^ OjiPL^ + ^iT^ D aPr^ ~ m.^ {'i> Pr + '^'^ Pl)^. (3.1) 

The 6-dimensional gamma matrices F"^ are defined in terms of the vielhein and of the flat-space 7° according to 
F"^ := e^7°. The covariant derivatives are defined as 

I'^A + lu^iWb - ^9Wa, Da ^Va + ^^fa,^, (3.2) 

^ y 

where uj'^ := |e°^Vy^eg are the spin connection with generators a^j^ := ^[73, 7b]- The symbol d acts on spinors 

according to "0 9 = ijjd(f> — {dil:)(j). Capital latin letters run from to 5 and denote the six dimensional space-time 
index while the hatted small Latin letters range from to 5 and corresponds to the fiat tangent six dimensional 
Minkowski space indices. 

We employ the left/right symmetric coupling scheme $ := r • u^Wa ~ 0; this form was extensively studied in 
(49I [soj and was applied to the six dimensional brane physics with warped geometry [l^ . 
The vielbein is defined through gAB = e^eas = Vah^A^-^B- We use the following form 

e°^M{p), ei^M{p)e"»'5i, fc = l,---,3, 
Cp — COS0, ~ sinfl, 

ei^-l{p)sme, e^g^l{p) cose . (3.3) 
A straightforward calculation shows that the nonvanishing components of the corresponding spin connections are 

uj°^ = -cos9M', cjf = -cos6iAf'e-f^°*(5^, fc = l,---,3 

, ,05 _ „;„/3)\/r' , ,a5 _ OAT' „Hat xa 
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sinflAf, = - sin6lA/'e/"«*(5g, fc = 1, • • • ,3 
1 - I'. (3.4) 



The standard (Dirac-Pauli) representation of the gamma matrices in six dimension is given by 

r 11 \ .A._( -Hi \ .,5 _ f -h 



where 7'' are the standard representation of the usual gamma matrices. The six dimensional spinor can be 
decomposed into a four dimensional and an cxtradimensional components as ^'(x^, p, 6) = il;{x'^)(E)(Ui{p, 6), U2{p, 9))^ , 
where ip{x^), Ui{p,9) are four, two components spinor, respectively. In terms of the decomposition, the equations for 
the extradimensional components are 



-Me" 



{dp \de 2^ - U2 + MruT ■ nU, = wU^ 



Me^' (dp + \de - 2^ - ^^^) U, - Mmr ■ nU^ = WU2 (3.6) 



where 'ip{x^) is the solution of the 4 dimensional Dirac equation 
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Therefore, the Dirac equation in six dimensions reduces to a two dimensional eigenproblem where the eigenvalue w is 
the masses of the fermions measured on the four dimensional brane. The Hubble constant works as a time component 
of the vector potential vq '■— fVo which simply shifts the energy as E ^ E + vq. In terms of the change of spinor 
U = (f/i, into u = (mi,W2)^, 

/ 1 1 r dr' \ 

w(r,0) :=cxp (2 In A/(r) + - In L(r)--y -^^U{r,0), (3.8) 



the eigenproblem then becomes 



Hu = wu, (3-9) 



where the hamiltonian is given by 



rfiT ■ n — e ^^{dr—i—de, 
\dr + j — mr • n 

Note that we have introduced dimensionless coupling constant and eigenvalue 

m := 1/ — m, w := ^ — w (-J-llj 

\ K2 \ K2 

and also L := a/ K2/K4/, see (|2.7p . 

The hamiltonian p.lOp is invariant under time-reversal transformation defined by T := 27^ (E) t'^C, where C is 
the charge conjugation operator. One can also easily confirm that the hamiltonian H commutes with "grandspin" 
operator given by 

Ks-.^k + ^+n^, (3.12) 

where I3 := —i-^ is the orbital angular momentum in the extra space and where we introduced 7^ := I <S) for 
convenience (TB]. As a consequence the eigenstates are specified by the magnitude of the grandspin, i.e., 

= 0,±1,±2,±3--- , for odd n 
13 5 

/•iTg = ±-,±-,±-, • • • , for even n. (3.13) 
and it follows from the time-reversal symmetry that one finds that the states of ±1^3 are degenerate in energy. 



B. Asymptotics 



The general form of solutions to 



is given by 



/ (7i(r)e*(^'^3-5-f )« \ 
c,2(r)e*(^^"5+t)e 
/ii(r)e*(^3+i-f)e 



(3.14) 



Far from the origin, the functions follow an exponential behavior that depends on /3; for /3 = 0, 



9i 

92 



^e"^'', hi ^ e~ 



(3.15) 
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where 



A 

V -- 

B 



1 



1 



2Lq Lq 
1 

2rhLi 
1^1 

2Lo Lq 

1 



rhLo 



2min 



filLr] 



K3 + 



rh^Ll- — 



K 



^-3 i- 7: - TT 



3+0 + 77 



For /3 < 0, the solutions arc simply gi,hi ^ e 

Close to the vicinity of the origin, the regular solutions to the linearized equations are 



92 



(3.16) 



(3.17) 



Therefore, normalized solutions of the eigenequations 
exist. 



interpolating between the near origin and far region should 



C. Schrodinger type equation 

The eigeneqaution p.9p can be recast into a set of Schrodinger-like second order differential equation. If we eliminate 
the components (hi, /12) , after a lengthy calculation we finally get 



- (Mf + ^)a, + ii + ^„ Widr + Vi2 \ ( Ui 



where I := i^s — \ — V^. Here we used the following replacement 
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32/ y g-| J"(A22 + i-fl)dr'y^ 

Similarly for the lower components ft-i, /12, we get 

/^-^r-(Tf + x)^'- + ^+t^ll ~W^29. + C/l2 \ /I'l 

1^ -w^ia. + t/2i -a2-(Mf + ^)a, + ^^±^ + c/22; V2 

The explicit forms of Aij, Uij, Vij, Wi] i,j = 1,2 are summerized in Appendix A. 



(3.18) 



(3.19) 



(3.20) 



D. Spectral flow 

Instead of solving the second order differential equations p.l8p . p.20p . we directly find solutions of the original form 
of the Dirac equation p.9p . The computational method is based on a plane wave expansion of the spinor and the 
matrix diagonalization scheme. (The detail was described in Sec. Ill of our previous paper p^.) 

Interesting property of such isolated bound states is that they dives from positive energy to negative if background 
fields change. This is called the spectral flow or the level crossing picture [i^. The spectral flow is defined as the 
number of eigenvalues of Dirac Hamiltonian that cross zero from below minus the number of eigenvalues that cross 
zero from above when varying the properties of the background fields. According to the index theorem, a nonzero 
topological charge implies zero modes of the Dirac operator [46l | . The number of flow coincides with the topological 
charge and zero modes emerge when the flow crosses zero. The level crossing picture was extensively studied in 



the Dirac equation with non-linear chiral background l45l. ISll. with Higgs field in the Abelian-Higgs model [52|, 155, 
and with non-trivial gauge fields fe.g.jnstanton.meron) |54l456l |. The mechanism can be interpreted as a quantum 
mechanical description of fermion creation/annihilation. 

In topological charge = 3, the index theorem states that three positive energy levels should dive to negative 
continuum, as shown in Fig. [S]for (to, n) = (1, 3) as well as for (to, n) — (3, 3). However, in the latter case, the picture 
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FIG. 5: The spectral flow with (m, n) = (1, 3),a = 0.1, 7 = 0.002 and p = 0.0 (left), and (m, n) = (3, 3), q = 0.001, /3 = 0.0, 7 = 
0.00125 (right). For (m,n) = (3,3), the bold letters beside the lines indicates number of the degeneracy. 



is much more sophisticated than for m = 1. A doubly degenerated heavy state interchanges with a light single state 
at junction A. The junction B is more complicated. At least, two isolated levels and two doubly degenerated levels 
interact each other. (Similar behavior has already been observed in a somewhat different context in [52|.) After the 
spectral flow, one easily confirms that three levels dive from positive continuum to the negative. In the following 
analysis, we concentrate on the simpler case (m,n) = (1,3) because it is more easy to get physical intuition. Figl6] 
shows the localization properties of components of the spinors and the scalar density which is defined as 

2 

^oW = II / deu\{r,e)u,{r,e) (3.21) 

i=l 

for several coupling constant corresponding to {m,n) = (1,3). We could produce the infinite tower of the spectra 
corresponding to K3 = 0, ±1, ±2 • • • as well as the complete set of the wave vector. Many of them form continuum 
(conducting level) and the lowest a few states are occupied levels. For small coupling constant, no localized solution 
exists. Increasing the coupling constant, the first degenerate states begin to localize and a single state follows. 
However, for a sufficiently large coupling constant, the spectra merges to the negative continuum and localized modes 
disappear. 



IV. HUBBLE PARAMETER AND THE MASS DIFFERENCE OF THE QUARKS 

The value of the Hubble parameter at the present epoch is about [13, [sH] 

Ho - 70 km s-^Mpc-i = 1.5 x 10"^^ MeV (4.1) 

In our previous analysis p^ . we used Skyrme parameters K2, ^4 of order y/ K2/K4 ~ 10^ MeV, so the dimensionless 
parameter 70 is around 

^0 = 1±H^ ^ 10-8 X (1.5 X 10-39)2 = 2.3 X 10-86, (4.2) 

which is apparently negligible. However, during inflation the Hubble parameter should be quite different, i.e., the 
typical value is around TJinf ~ 10^°GeV, which corresponds to 7inf = 10^* [59| . 

In this section, we discuss the behavior of the massive fermion levels (quarks) for changing 7. Since we have 
several model parameters, we need to fix some of them in an appropriate way. We always get the localized modes for 
first/second degenerate levels (which corresponds to states K3 — ±1) by suitably adjusting the coupling constant m 
for each 7. If we simply set zero for these levels, the third generation (the K3 = state) also becomes localized mode. 
So we compute the mass of the third generation as the level of this localized mode Cjks^o- In Fig 17] we show the result 
for the case of /i = 0.0, 0.1, 0.5. The mass differences monotonically increase as Hubble parameter 7 grows. Note 
that the solutions with n = is naturally interpreted as corresponding to the case oi kq = 0, which is excluded at 



9 




FIG. 6: The spinors (left) and the scalar densities (right) for the skyrmion backgrounds with (m,n) = (1,3), a = 0.01, /3 = 
—0.01,7 = 0.005 and ji = 0.1 for the several coupling constant m. A few of the components of the spinor sometimes are 
degenerate. 
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FIG. 7: The mass difTerence with n = 3, a = 0.05, /3 = 0. 
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FIG. 8: The spectral flow with n = 3, a = 0.0, /3 = 0.0, /i = 0.1, as a function of the cosmological constant. 
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FIG. 9: The effective potential V22, Uu for several parameters 7 with the case of n = 3, a = 0.0, /3 — 0.0, n = 0.1. 



least in flat space by Derrick theorem. Here, we flnd solutions for /3 = 0, it follows that in this case /i = is to be 
regarded as other limits, compaptible with Derrick theorem. Such a limit can be achieved for ^2 — >■ 00 or K4 = 0. 

Roughly speaking, smaller values of 7 corresponds to late stages in the evolution of our Universe, in which case the 
mass difference monotonically decreases. In this sense, the parameter 7 allows to follow the properties of the fermions 
at different epochs during the inflation. However, our previous analysis (loj showed that there exists a maximal value 
of 7 where we could find solutions. Typically the maximal values was of order 1 which is extremely far from 7inf . It 
follows that there are two ways of interpreting our model: 

(i) our brane solution is suitable for describing only of the late stage in the evolution of our Universe; 

(ii) our brane solution could describes all age of Universe, but different values of the parameters ^2,^4 should be 
employed for each epoch. 

The first possibility (i) is the most natural interpretation. The value of 7crit indicate that the solution starts with 
a late stage, which somehow makes sense since the model should be an effective model of some more fundamental 
mechanism. Therefore the baby-skyrmion can only describes the quite recent Universe, already reducing speed of the 
inflation. Gap of the fermions between first/second and third generations gradually decreases as time grows (FiglTl). 
The change is, however, quite subtle. It depends on /i, but the difference at the beginning is at most ten times larger 
than now. 

The second story (ii) seems strange. It is based upon a key assumption: the model parameters of the baby- 
skyrmion are function of the age of the Universe. We suppose that at early (or the beginning) stages, the value of 
7i,if = 7crit = 2.8. Therefore we estimate the model parameters at this epoch as 

= 2.8 = ^ X 10^6, ,. ^=2.8x10-26. (4.3) 

K2 K2 

In order to understand the structure of the brane in early time, it may be helpful to use known asymptotic solution 
of the baby-skyrmions in flat space of the form [TJ] '■ 

f[r) ^ '^%-i/2e-yw^ ^ (4.4) 



^From (|2.8p and one easily see 



= I!^/, = 3.6 X lO^V (4.5) 

K2 K4 



thus f ^ for this limit. This clearly indicates that the brane shrinks for the case of flnite fi, i.e. at early stages of 
inflation. The dimcnsionful mass difference uj becomes 

'"^w ^ 3.6 X lO^^w (4.6) 
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which indicates that the third generation blows up. It rapidly decreases as time grows (K2/K4 decreases) and finally 
arrive to a finite value at our epoch (1/^2/^4 ^ 10^). 

In conclusion, the mass difference decreases as time marches on, but in the picture (i), the effect is moderate so that 
the third generation would be observable even at the early stage while in the picture (ii) it would never be observed 
because the energy is almost over ^ 10^° than scope of any experimental facility. More detailed analysis is however, 
outside the scope of our present model. 

Next we investigate the broader range of the parameter 7. If we choose the coupling constant rh = 0.62 for a 
solution of a = 0.0, /? = 0.0,^ = 0.1, it makes the first /second levels to zero at 7 = 0. We could observe different 
type of the spectral flow (Fig|S]). For decreasing value of 7, the levels = 0, ±1 dive from positive continuum and 
become localized modes. In this parameter case, the level = ±1 begin to localize to the origin at 7 ^ 0.09 and 
when 7 ^ 0.03, the level = becomes the localized mode. In FiglH we plot diagonal part of the effective potentials 
V22, Uii (the explicit forms are shown in the appendix, (|A4|) . (|A7p .) Those exhibit the volcano form. For larger value 
of 7, the well becomes deep in depth but narrow in size. At a critical point (e.g. ,7 0.09), the potential suddenly 
blows up and after that, no any localization modes appear. 

This mechanism explains creation of the particles (the massive fermions) in an expanding Universe. Generally, the 
particles appear by the pair creation/annihilation process. In the level crossing picture, the particle pair creation 
occur during the level crossing the zero from negative to positive. In our case, however, the expanding brane captures 
and localizes the fermions in the extra dimensional space time, which means that during inflation massive fermions 
suddenly appear on the brane (from upper continuum) and become lighter as 7 reduce, i.e. as inflation continues. 



V. MAXWELL FIELD ON AND OFF THE BRANE 



In this section, we consider problem of a vector field in the background of the brane. According to the analysis 
for the 2+1 baby-Skyrme model with a U(l) gauge field [g^l, we introduce the action with the gauge field in six 
dimensions 



^Maxwc! 



n = / d^x 



-9 



4e 



MN 



(5.1) 



and the usual derivative in the baby-Skyrme action 5'baby is replaced to the covariant derivative T^m = Om + 
AmUoo X . with the U(l) gauge field Am- The Tmn '■= QmAn — OnAm is the Faraday tensor in six dimensions. 
Variation of the action S :~ 5'baby + •S'Maxwcii with respect to A]\j leads to 



y Mr = eJ , 



= {K2(noo X n) • (V^n) - K4 [(rioo x n) x VMn] . [V^n x V'^'n] } 
We employ following paramctrization for the gauge field 

Am ^ iA^{xnF{r,e), A.,{r,9)), 



(5.2) 
(5.3) 



where ^, v refers to the four dimensional indices, i.e. = t, x, y, z and i means the extradimensional components 
such as i,j = r,0. Using this parametrization, the Faraday tensor can be expressed as 



f_ii/ F 


~A^d,F 


A^d,F 





(5.4) 



where F^^) ~ dfj,Au{x) — duA^(x) is the four dimensional Faraday tensor. 

Note that the action (|5.ip contains a term jgF which clearly means that not only the Planck mass, 

but also the electromagnetic coupling acquire the hierarchy. The effective four dimensional electromagnetic coupling 
would then become 



64 



/F2(r, e)L{r)drde 



(5.5) 



In fact, as we shall see in the following, the equation for F turns out to be linear, thus we can choose the normalization 
of F such that / F{r,9 fL{r)drde = 1. 

The equations for the gauge fields in the background of the baby skyrmion then become 



-(V'^A^(a;))V,i? + V^J-jf =eJ„ 



(5.6) 
(5.7) 
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The four dimensional components of the source is actually proportional to say = J(r, 9)A^ where J can 

be estimated via the extradimcnsional components, i.e. /, /', M and F. The equations for the extradimensional 
components then reduce to 

DF ^ eJ ^ H^F, (5.8) 

S7'Fj, - eJ, = \d,F, (5.9) 

where we used 

V^J-W ^ -f,^A,, (5.10) 

V'^A^(x) = A, (5.11) 

where A, /i are the variable separation constants. Quite interestingly, as a consequence of the variable separation (|5.3p . 
above equations naturally implement the Maxwell equations and a gauge condition. ^From a four dimensional point of 
view, the vector field should be massless, we impose /i = and then (|5.10p is the four dimensional Maxwell equations. 
Also if we choose a special gauge choice A = 0, (|5.1ip is exactly the Lorentz gauge condition in four dimensions. 

Note that direct computation shows that depends on the four dimensional components only through a pro- 
portionality factor with and Ji depends on the extradimensional coordinates only. Note also that the covariant 
derivative acting on four dimensional objects reduces to the four dimensional covariant derivative. 

We can successfully decouple the four dimensional and extradimensional part when we impose the four dimensional 
Maxwell equations and the Lorentz gauge condition. We parameterize the extradimensional part as Ai{r,9)dx'' = 
nAg(r)d6 and the extradimensional dependence of the four dimensional components as F{r,9) = F{r); which lead to 
the following two decoupled equations 

0, (5.12) 
- e + 1) sinV = 0. (5.13) 

For the boundary conditions we impose 

F'{0) = 0, Ag{0) = 0, F{0) = 1, Ae{oo) = 0, (5.14) 

where the first two are the regularity condition of the equations, the third one is an arbitrary choice of normalization 
and the last should be imposed in order to have no flux at inflnity. 

We should stress that all parameters do not lead to localizing gauge fields. The condition for localizing mode is of 
course the function F decays to zero as r grows. Since we impose that the four dimensional components are massless, 
the only parameter we can vary is actually the electromagnetic coupling. 

For solving the equations, we use a standard fourth order Runge-Kutta method, with shooting for (|5.12p and the 
backward integration for (|5.13p . The reason for using backward integration is that it avoids the need of shooting to 
find the proper decay. As we shall discuss later, the source term in (|5.13p is non-trivial, and then it imposes the value 
of the function Ae at the origin. 



F" 
A'' 



2^ 


V 


M 


+ L 




L' 


M 





F' + eF + 1) sin^ / = 
A', + eAe (/'^ + l) sin^ / - 



A. Four dimensional gauge field conflnement 

The near origin expansion of the four dimensional form function F is given by 

4(?i + ly 

where Fq is a real constant and /„ is the coefficient of the r" term in the near origin expansion of the function /. 

We focus on the case of vanishing bulk cosmological constant. In this case, for large values of r the function F 
behaves like 

F{r)=Fl + ^ + o{^^\ (5.16) 
where F^^, F^ are real constants that depend on the parameters, especially on e. 
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a = .05, ^ = 0, fi = 0.1, m = 1, n = 3 



e 




FIG. 10: The values of e as a function of 7 for the first three locahzing modes. 



a = 0.05, J3 = 0,y = 0.001, /i = 0.1, m = 1, n=3, e = .019867 a = 0.05, /3 = 0. y =0.001. /i =0.1. m = 1, n=3, e = .976949 




FIG. 11: The extradimensional dependance of the four dimensional components of the gauge field. We added the best fit of 
the decay, starting from the value r = 15. 

In principle, there should be a massive tower of four dimensional gauge fields. However, the lowest state is expected 
to be massless, so we concentrate on the zero mode. /^From the fact is that not all values of the coupling lead to 
localized zero modes, i.e. modes for which the function F — > for large values of r (in other words for which = 0). 

Our result indicates that only specific values of e lead to localized four dimensional gauge fields. The process is 
somewhat similar to that of an eigenvalue problem, and in this case the coupling constant plays a role of the eigenvalue. 
In fact, when looking at (|5.12p . up to the sine terms, the equation certainly looks like an eigenvalue equation. The 
fact provides a clever quantization mechanism for the elementary electric charge due to the extradimensions. 

Furthermore, we find that values of the electromagnetic coupling for the localized modes depend on the inflationary 
parameter. As discussed in the previous sections, varying the Hubble parameter in the model can be interpreted as 
looking to different time slices of the universe. Seen in this interpretation, our results suggest that the value of the 
electric charge can evolve with time. This is illustrated on Fig. [101 We also show two particular localized modes for 
specific values of the parameters in Fig. 1111 



B. Gauge field in the extradimensions 

In this subsection, we consider the extradimensional component of the gauge field Ag in (|5.13p . The main difference 
of the equations is (|5.13p contains a source term. As a consequence, we cannot invoke the linearity to choose 
normalization, instead, the source term imposes the value at the origin. One should keep in mind that we study the 
case of gauge field in the background of the baby-Skyrme. As a consequence, the results presented here might be 
affected by backreactions, especially in the case of larger electromagnetic couplings. 

The near origin behaviour of the function Ag is given by 

Ae^ao + a,r' + + O {rr+' , (5.17) 

n(n + 2 
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a = 0.05, /3 = 0,y = 0.001, = 0.1, m = 1, n = 3 




FIG. 12: The value at the origin of the extradimensional component of the gauge field as a function of the electromagnetic 
coupling. 
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FIG. 13: The profile of the extradimensional component of the gauge field for three specific values of the coupling and specific 
values of the background parameter. The larger the coupling, the more the field is delocalized. 

where ao,a2 are real constants and fn is again the coefficient of the r" term in the near origin expansion of the 
function /. In the case of vanishing bulk cosmological constant and non-vanishing brane cosmological constant, the 
field Ae decays according to 

As = ^+o(^y. (5.18) 

Wc find that for small values of the electromagnetic coupling, the extradimensional gauge field is more localized 
around the brane. As the coupling increases, the gauge field starts to produce a finite number of oscillations around 
some positive values before decaying to zero. We note that the number of oscillations increases with the electromagnetic 
coupling. Furthermore, the more the field oscillates, the more it is delocalized (i.e. the size grows). Note that the 
gauge field does not develop modes as it was doing in the four dimensional case; it always stays positive. We illustrate 
the above discussion on Figs. [T2]and[T3l 
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VI. SUMMARY 



We have studied fermion localization on the baby-skyrmion brane with positive cosmological constant. We especially 
have concentrated on the solutions with topological charge iV = 3. The matter distributions have a ring shaped but 
for m = 3 it is more concentrate on a region close to the origin. The fermions coupled with the baby-Skyrme field 
has an isospin doublet structure, thus we naturally could take into account the flavor symmetry of the quarks. In 
terms of the index theorem, number of the fermions localized on the brane coincides with the topological charge 
N. For TV = 3, we obtained the three localized modes during a parameter range of the coupling constant to. The 
doubly degenerate levels are always lower than the single excited state, the levels are occupied by quarks/leptons of 
lower two generations. Especially we studied about the mass difference between the first two generations and the 
third generation for changing dimensionless Hubble parameter 7. The behavior exhibits typical spectral flow. For 
decreasing value of 7, the three localized modes appear from the upper continuum when 7 ^ 0.03, which indicates 
new particle creation mechanism. 

Our scheme is symmetric to the fiavor (isospin) degrees of freedom thus (w, c, t) sector and {d,s,b) sector are 
degenerate. In order to split the degeneracy, we introduce explicit symmetry breaking of the skyrmion parameters 
K2, K4. For (d, s, b) sector the choice K2/K4 ~ 10^ MeV reproduce the mass difference between the first/second and 
the third generationsm while for (u, c, t) sector we adopt -y/ K2/K4 ~ 10^ MeV to achieve an agreement with observation. 
This prescription equals to introduce the flavor hierarchical structure to the coupling constant m (for a fixed rh) in 
terms of the relation p. lip . On the other hand, in Ref.flBl one of the authors of the present paper considered 
the breaking of the time-reversal symmetry by taking into account the rotational symmetry breaking of the baby- 
skyrmions. The authors of [l6j introduced an arbitrary deformation parameter by hand which successfully reproduce 
the mass splitting of the first /seco nd g enerations. Recently we have found the existence of the Z2 symmetric, deformed 
brane in the baby-Skyrme model [6l[. It is worth to computing the spectra of the fermions coupled to the deformed 
brane. As a conclusion, for the flavor hierarchy we could solve by introducing the hierarchical coupling constant of the 
fermions and the baby-skyrmions. Also the geometrical structure of the brane (background baby-skyrmions) brings 
about a solution to the mass hierarchy within the family. A more quantitative analysis implementing all these issues 
will be discussed in the forthcoming paper. 

In this paper, we employed left/right symmetric coupling scheme of the fermions and the baby-skyrmions. And we 
used the six dimensional generalization of the Dirac-Pauli representation for the gamma matrices. Although there are 
several advantages to computing the massive modes, the mechanism of localization of chiral fermions on the brane 
is still absent. The study is mandatory to achieve full understanding on properties of our SM particles and will be 
reported in near future. 

On another hand, we studied the localization properties of gauge fields and found that indeed, it is possible to 
localize four dimensinal gauge fields on the brane for very specific values of the electromagnetic coupling. We studied 
the gauge field in the background of the brane, i.e. without backrcactions. Although we believe that it is already 
quite instructive for the four dimensional part of the vector field, we think that the cxtradimensional components of 
the gauge fields should be considered with full backreactions. The good news is that, at least in our approach, the 
cxtradimensional sector and the four dimensional sector of the vector field decouple. 

We forced massless gauge fields, even if it is well known that higher dimensions usually have the effect of producing 
mass towers for the fields of interest. In our case, we know that electromagnetic gauge fields are massless from a four 
dimensional point of view. This implies that we should at least have a massless mode for the four dimensional gauge 
fields. This zero mode should actually be the first state of a massive higher energy tower of vector fields (which we 
do not study in this paper). 

As an plausible extension of the present study is to consider non-abelian gauge group confinement on solitonic branes. 
As we observed in this paper, imposing existence of massless vector mode automatically introduces quantization on 
the coupling, which is somehow the inverse of the usual approach, namely, compute masses for a given coupling in 
order to find localized modes. This is certainly valid for a low energy approximation of massless vectors, such as 
U{1) fields. However, for vector bosons like SU(2), provided the W, Z particles, it does not work because they have 
observed masses. One should then either recall the standard procedure of computing the mass for given coupling or 
find a new mechanism that generates mass and predicts coupling. 
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Appendix A: The components for the effective potential 



The Vij,Wi = 1, 2) in (|3.18p are summarized as 

Vii = -A' + -Mil + - + All- - —r 



2^1-2— -2——--— + -—--— ' 



M M2 ML 2 L 4 2 £2 ' 7,2 



^--^-'-21^-+lJ+^/ ' 2L L 



1^22 = x^',, + TU22 + ^ +A22— 



2 



™2 



2? = —^oi n ^22 n + /I22 

22 ^ 22^ ^\ 22^ ^ 22 ^ 

~ ~m"" 17 17 " 2~L ^ 4"l2" " 2I2 Z2 

For p.20p . the components t/^ , (i, j = 1,2) are written as 

^"-2^22 + 4(^22 + ^) +^22^-^+m^ 

_2^_2— ~2— — -i— G + l)(L^-l) 

~ ~M ~ M2"~ a7"l " 2I7 ^ 4L2" " 2Z2 + J2 ' 



. , M' IL' l + n+l\ in A a w.' 

U12 = Ai2\ Uii + - +2 + + 6-2.1 (^ii-^22)dr 

2\ l) M 2 L L 

C/21 = A21 ( 4(A22 + i)+2^ + l^ + i±l]e-'^riA..-M.)'ir' 



2\ LI M 2 L 

[/22 = -A'li + -Mil + — +^11 ^+m2 

_2^-2— -2— — -i— ----- (^ + ^t + l)(-^'- 1) 

~ ~m" ~ M2" " ^ 2ir ^ 4 L2" ^ 2I2 £2 

The Aij{i,j — 1,2) are obtained as 

/ a;2 2\~^/ 2" 2 il/'a;2 ojto M'ujrh \ 

/ ^2 „\-^/mu} fn M'\ fnfii ^f'\ 2 / 



A21 ■= - — ^ - m sin f \ m cos f sin f m f' 

VM2 J \AI ^\L m) ■' \ L ■' M ' ■' 

. / ~ 2\ V ~ 2 • 2 /• , M'u? ujifi . , M'ujrh > 
^22 — T - m ] m —sm f H 5 sm f f H ^ cos f , 

and the Ai2,A2i can be estimated similar to Aij, and are found to be 

A12 := [jj-2 m^) sm/(- + — ) - meos/(— sm / - — ) - f 

^21 := - — ^ - sm f + m cos f sm f — - f 



18 



[1] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) 'ar Xiv:hep-ph/9803315, 

[2] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Rev. D 59, 086004 (1999) a rXiv:hep-ph/9807344| 

[3] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) arXiv:hep-ph/9905221 

[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) arXiv:hep-th/9906064 

[5] A. G. Cohen and D. B. Kaplan, Phys. Lett. B 470, 52 (1999) arXiv:he p-th/9910132| ' 

[6] R. Gregory, Phys. Rev. Lett. 84, 2564 (2000) arXiv:hep-th/9911015 



[7] T. Gherghetta and M. E. Shaposhnikov, Phys. Rev. Lett. 85, 240 (2000) [arXiv:hep-th/0004014 

[8] M. Giovannini, H. Meyer and M. E. Shaposhnikov, Nucl. Phys. B 619, 615 (2001) 'arXivihep-th'/OlOiTTS 



[9] C. Ringeval, P. Peter and J. P. Uzan, Phys. Rev. D 71, 104018 (2005) arXiv:hep-th/0301172 



[10] E. Roessl and M. Shaposhnikov, Phys. Rev. D 66, 084008 (2002) [arXiv:hep-th/0205320| . 

[11] T. H. R. Skyrme, Proc. Roy. Soc. Lond. A 260 (1961) 127. 

[12] L. D. Faddeev and A. J. Niemi, Nature 387, 58 (1997) (arXiv:hep-th/961 0193 . 

[13] B. M. A. Piette, W. J. Zakrzewski, H. J. W. Mueller-Kirsten and D. H. Tchrakia n, Phys. L ett. B 320, 294 (1994). 

[14] B. M. A. Piette, B. J. Schroers and W. J. Zakrzewski, Z. Phys. C 65, 165 (1995) [arXiv:he p-th/9406160 . 

[15] A. E. Kudryavtsev, B. Piette and W. J. Zakrzewski, Eur. Phys. J. C 1, 333 (1998) 'arXiv:hep-th/961 1217| . 

[16] Y. Kodama, K. Kokubu and N. Sawado, Phys. Rev. D 79, 065024 (2009) ;arXiv:0812.2638 [hcp-th]]. 

[17] Y. Brihaye, T. Delsate and B. Hartmann, Phys. Rev. D 74, 044015 (2006) arXiv:hcp-th/0602172 . 

[18] Y. Brihaye and T. Delsate, Class. Quant. Grav. 24, 1279-1292 (2007) arXiv:gr-qc /0605039_ 

[19] Y. Brihaye, T. Delsate, Y. Kodama, N. Sawado, Phys. Rev. D82, 106002 (20 10r^arXiv:1067.0736l [hep-th]]. 



[20] A. Kehagias and K. Tamvakis, Phys. Lett. B 504, 38 (2001) arXiv:hep-th/0010112 



[21] C. Ringeval, P. Peter and J. P. Uzan, Phys. Rev. D 65, 044016 (2002) arXiv: hep-th/0109194l. 

[22] A. Melfo, N. Pantoja and J. D. Tempo, Phys. Rev. D 73, 044033 (2006) arXiv:hep-th/0601161|. 

[23] R. Koley, J. Mitra and S. SenGupta, Phys. Rev. D 78, 045005 (2008) arXiv:0804.1019 [hep-th]]. 

[24] Y. Hosotani, S. Noda, Y. Sakamura and S. Shimasaki, Phys. Rev. D 73, 096006 (2006) arXiv:hep-ph/060 1241| 

[25] Y. X. Liu, H. T. Li, Z. H. Zhao, J. X. Li and J. R. Ren, JHEP 0910, 091 (2009) arXiv:0909.2312 [hep-th]J 

[26] H. Guo, A. Herrera-Aguilar, Y. X. Liu, D. Malagon-Morejon and R. R. Mora-Luna, arXiv:1103.2430 [hep-th]. 

[27] K. Agashe, A. Falkowski, I. Low and G. Servant, JHEP 0804, 027 (2008) arXiv:0712.2455 [hep-ph]]. 

[28] S. Randjbar-Daemi and M. E. Shaposhnikov, Phys. Lett. B 492, 361 (2000) arXiv:hcp-th/ 0008079] . 

[29] M. V. Libanov and S. V. Troitsky, Nucl. Phys. B 599, 319 (2001) arXiv:hep-ph/0011095j . 

[30] A. Neronov, Phys. Rev. D 65, 044004 (2002) arXiv:gr-qc/0106092 . 

[31] S. Randjbar-Daemi and M. Shaposhnikov, JHEP 0304, 016 (2003) arXiv:hep-th/0303247]. 

[32] S. Aguilar and D. Singleton, Phys. Rev. D 73, 085007 (2006) arXiv:hep-th/0 602218^. 

[33] M. Gogberashvih, P. Midodashvih and D. Singleton, JHEP 0708, 033 (2007) larXiv:0706.0676l [hep-th]] 



[34] S. L. Parameswaran, S. Randjbar-Daemi and A. Salvio, Nucl. Phys. B 76 7, 54 (2007) larXiv:hep-th/ 0608074 

[35] L. Zhao, Y. X. Liu and Y. S. Duan, Mod. Phys. Lett. A 23, 1129 (2008) larXiv :0709.r520 [hep-th]]. 

[36] Z. q. Guo and B. Q. Ma, JHEP 0808, 065 (2008) arXiv:0808.2136 [hep-ph]]. 

[37] Z. Q. Guo and B. Q. Ma, JHEP 0909, 091 (2009) arXiv:0909. 4355 [hep-ph]]. 

[38] V. A. Rubakov and M. E. Shaposhnikov, Phys . Lett. B 125, 139 (198 3). 

[39] K. Akama, Lect. Notes Phys. 176, 267 (1982) [arXiv:hep-th/0001113| . 

[40] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147, 277 (1979). 



[41] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D 61, 033005 (2000) arXiv:hep-ph/99 03417 

[42] G. R. Dvah and M. A. Shifman, Phys. Lett. B 475, 295 (2000) arXiv:hep-ph/0001072 . 

[43] P. Q. Hung and M. Seco, Nucl. Phys. B 653, 123 (2003) arXiv:hep-ph/0111013 . 

[44] J. M. Frere, M. Libanov and F. S. Ling, JHEP 1009, 081 (2010) arXiv:1006.5196 [hep-ph]]. 

[45] S. Kahana and G. Ripka, Nucl. Phys. A 429 (1984) 462. 

[46] M. F. Atiyah, V. K. Patodi and I. M. Singer, Math. Proc. Cambridge Phil. Soc. 79, 71 (1976). 

[47] P. Eslami, W. J. Zakrzewski and M. Sarbishaei, arXiv:hep-th/0001153 

J. C. U. Ascher and R. D. RusseU Math, of Comp. 33, 659 (1979) 

[49] T. Jaroszewicz, Phys. Lett. B 146, 337 (1984). 

[50] M. S. Carena, S. Chaudhuri and C. E. M. Wagner, Phys. Rev. D 42, 2120 (1990). 

[51] N. Sawado and N. Shiiki, Nucl. Phys. A 739 89 (2004) , arXiv:hep- ph/0402084| . 

[52] Y. Burnier, Phys. Rev. D 74, 105013 (2006) [arXiv:hep-ph/0609028t 



[53] F. L. Bezrukov, Y. Burnier and M. Shaposhnikov, Phys. Rev. D 73, 045008 (2006) [arXiv:hep-th /05 1 2 143] . 

[54] N. H. Christ, Phys. Rev. D 21, 1591 (1980). 

[55] H. B. Nielsen and M. Ninomiya, Phys. Lett. B 130, 389 (1983). 

[56] J. E. Kiskis, Phys. Rev. D 18, 3690 (1978). 

[57] A. G. Riess, L. Macri, S. Casertano, H. Lampeitl, H. C. Ferguson, A. V. Filippenko, S. W. Jha, W. Li et al., Astrophys. 

J. 730, 119 (2011). [arXiv:1103.2976l [astro-ph.CO]] . 
[58] F. Beutler, C. Blake, M. CoUess, D. H. Jones, L. Staveley-Smith, L. Campbell, Q. Parker, W. Saunders et al., 

[arXiv:1106.3366 [astro-ph.CO]]. 
[59] I. Cho, Phys. Rev. D 69, 105019 (2004) [arXiv:hep-th/0402125| . 



19 



[60] J. Gladikowski, B. M. A. G. Piette and B. J. Schroers, Phys. Rev. D 53, 844-851 (1996) 
[61] T.Delsate, M.Hayasaka and N.Sawado, in preparation. 



20 



